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■ Abstract 
O , 

C , We consider the time-delayed coincidence counting of two photons emitted in a cascade by a 

o 

single particle (atom, molecule, nucleus, etc). The time-dependence of the probability amplitude 



(N 



O 

a: 



of the second photon in the cascade has a sharply rising leading edge due to the detection of the 
first photon, as results from causality. If a macroscopic ensemble of resonant two-level absorbers 

is. 

is placed in the path of the second photon between the radiation source and the detector, the 
photon absorption does not follow Beer's law due to the time-asymmetric shape of the photon. 



For very short delay times almost no absorption takes place, even in an optically dense medium. 



We analyze the propagation of such a second photon in a thick resonant three-level absorber if 
a narrow electromagnetically induced transparency (EIT) window is present at the center of the 
absorption line. It is shown that the EIT medium can change the asymmetric time dependence of 
the photon probability amplitude to a bell shape (EIT filtering). This bell-shaped photon interacts 
much more efficiently with an other ensemble of two-level absorbers chosen, for example, to store 
this photon and the information it carries. 

PACS numbers: 42.50. Gy, 42.50.Dv, 76.80.+y 
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I. INTRODUCTION 



In the last decade, there has been a considerable growth of interest in single photon 
experiments in the visible or near infrared domain, mostly related to information storage, 
quantum computing and quantum cryptography. An essential element in these experiments 
is a single photon source named "photon gun". It emits one and only one photon when 
the experimentalist "pulls the trigger". Faint laser pulses with extremely low mean photon 
numbers are only an approximation of single photon pulses. Their state is close to a Glauber 
state also containing two and more photons, and in the semiclassical description the laser 
pulse is modelled by a Gaussian pulse. Recently, several types of true single-photon sources 
were reported, which are based on the laser excitation of a single trapped ion or atom 111 , a 

q n 

single organic dye molecule in a solvent [2j , |3|] or a single nitrogen- vacancy color center in a 
diamond nanocrystal J] . A single-photon turnstile device that uses a single quantum dot or 
single quantum well was also proposed 15] . Such a photon is generated spontaneously from a 
single quantum object (atom, molecule or quantum dot) that is placed at time to = in an 
excited state by a short laser pulse or electronically by injecting a single electron and a single 
hole to annihilate in a light emitting domain (a central quantum well in a p-n junction). 

Single photon sources based on a single emitting particle can be divided in two kinds. A 
source of the first kind radiates in free space vacuum modes. For example, sources 
and |4| belong to the first kind. A source of the second kind radiates to a cavity mode. If a 
high-finesse cavity is used in a regime of strong coupling in cavity quantum electrodynamics, 
almost all radiation is collected in the active cavity mode and transferred through the cavity 
mirror (loss channel) in a well defined direction (see proposals based on the Jaynes- 

Commings model |8|). This process can be made deterministic and it allows the generation 

nnn 

of a single photon with a controlled waveform [9],[10J,[11J. In a quantum network consisting 
of spatially separated nodes connected by quantum communication channels it is preferable 
to operate with photons having a shape symmetric in time . In this paper we show that 
there is an additional physical argument demanding a time-symmetric photon wave packet 
if it is supposed to use such a photon in a atom-field interaction protocol for storage and 
retrieval of the radiation state. 

We consider a single photon source of the first kind, which intrinsically produces a photon 
with asymmetric temporal envelope because of causality. The resonant interaction of such a 
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photon with an ensemble of two- level absorbers (atoms, molecules, nuclei, etc.), referred to 
as a macroscopic absorber, is analyzed. We show that multiple scattering of the photon in 
the forward direction by the resonant particles displays unusual properties. They originate 
from the knowledge of the instant of time when the source particle was prepared in the 
excited state and then a photon was emitted, i.e., when the "trigger" of the photon gun 
was pulled. In the semiclassical language one can say that the photon wave packet has a 
particular time envelope having properties that are very unusual and different from those 
known for bell-shaped pulses. We compare the atom-field interaction of such a photon with 
that for a photon wave packet with a Gaussian time-envelope. Then, for a photon source of 
the first kind we analyze the possibility to reshape the photon field envelope with the help 
of electromagnetically induced transparency (EIT). The reshaped photon time envelope is 
close to the one with a Gaussian envelope and its interaction with another absorber (used, 
for example, for information storage or quantum computing) becomes the usual one typical 
for bell-shaped pulses. 

The paper is organized as follows. In Sec. II we discuss the influence of causality on 
the spectrum of the single photon emitted by a source of the first kind and show that it 
introduces a broad part. In Sec. Ill we consider the propagation of such a photon in a thick 
resonant absorber. In Sec. IV we study the photon filtering through a narrow EIT window. 
We show that the narrow spectrum part of the photon and the broad part are separated in 
time. Application of single photon filtering for quantum storage is discussed in Sec. V. The 
conditions for level mixing induced transparency for gamma radiation are considered in Sec. 



II. THE INFLUENCE OF CAUSALITY ON THE SPECTRUM OF A SINGLE 
PHOTON SOURCE OF THE FIRST KIND 

Assume that at time t = a single particle (generally a quantum object such as an 
atom, a molecule, a nucleus, etc.) is placed in an excited state e. If the radiation, emitted 
by this particle, has only natural lifetime broadening due to the coupling with the vacuum 
modes of free space, the radiation state is described by the usual expression [l^ 



VI. 
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where |lk) is a single-photon Fock state of the mode Uk with wave vector k, ro is the 
location of the particle, i>\. = u>k — u Q is the frequency difference of the fc-mode and the 
resonant transition from the excited state e to the ground state g, and 2A p h is the decay 
rate of the excited state e. The coupling parameter of the radiation with the source is g^. 
This state is actually a wave packet consisting of many Fock states and normalized such 
that in total it contains a single photon. Let 

(r, t) = jy^ a * e ~ Wkt+lkr > ( 2 ) 

k 

be the electric field operator containing only the annihilation operators a^, is the unit 
polarization vector, is a normalized amplitude of the mode k. Performing the sum over 
the wave vector k in the expression for the single photon field b(t) = (0| E^ + \r,t) \b), one 
obtains [3| 

b(t) = ^9 (t - d/c) e -(-o+A ph )(t-d/ c ) ; (3) 

Lb 

where d — |r — ro| is the distance from the source, So is a normalized amplitude and G(t) is 
the Heaviside step function. This wave packet has a sharply rising leading edge at t = d/c 
and an exponentially decaying tail. The former is defined by the time to = at which 
the source is placed in the excited state and the latter specifies the coherence time or the 
mean correlation time of the photon r p h = 1/A p h [14] . Such a time dependence of the 
single-photon field was detected from radiation of a single organic dye molecule in solvents 
j^j], in time delayed coincidence measurements (TDCM) of photons emitted in an atomic 



cascade (optical domain) [3| and in a nuclear cascade (gamma domain) ll£ , 12 , 3- Q • 
A typical energy diagram and a detection scheme of TDCM are shown in Fig. l(a,b). 
The delayed coincidence counting technique leads to the detailed observation of the time 
correlation between successively emitted photons a(t) and b(t) in the h — > e —>■ g cascade. 
The photons of interest, a(t) and b(t), can be selected by interference filters (optical domain) 
jl^ l or these photons can be of very different ener gies such that different kind of detectors are 
used to detect them separately (gamma domain) |l^.llfil l^IQ, 3]- In the TDCM detection 
scheme the energy of these photons is converted by high-speed, high-gain photomultipliers to 
electric pulses, which are then amplified and sent to fast discriminators. A time-to-amplitude 
convertor is turned on by a pulse due to a a(t) photon and shut off by one due to a b(t) 
photon. 
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Actually it is the probability P(t) = \b(t)\ 2 and not the complex probability amplitude 
b(t) which is measured. When a photon is detected at time t = t n , we have a collapse of the 
photon wave function at the detector place, which shows as a single count of the detector. 
After many events (detector counts), one can reconstruct the actual time distribution of 
the probability P(t). One can also extract information about the absolute value of the 
probability amplitude since \b(t)\ = \JP(t). If we know that at time to = the source 
particle is placed in the excited state e and there is no absorber between the source and 
the detector, this probability is zero for t < 0, it has a maximum, 1, at t — +0 and from 
then on it decreases exponentially. The time resolution of this measurement is defined by 
the response time of the detector and electronic circuit. Therefore, the stepwise rise of the 
probability P(t) at t = is usually slightly smoothened. 

For many applications it is important to know how a single photon field interacts with 
either a single absorber or an ensemble of quantum absorbers. In the time domain the 
propagation of such an asymmetric field through a thick absorptive medium was considered 
quantum mechanically and semiclassically for gamma photons and for small ampli- 
tude pulses in quantum optics [21]. It was shown that if the source and the absorber have 
the same resonant frequency and the same linewidth, the transmitted probability amplitude 
is b(T,r) = exp(—iu T)b (T,T), where 

6o(T, r) = e~ A ^J (2 ^TK~^) 0(r), (4) 

is the time envelope normalized to 1 for T = and r = 0, r = t — l/c is the local time, I is 
the physical length of the macroscopic absorber, T = aol/j is its effective thickness, Jq(x) is 
the zero-order Bessel function, and 0(r) is the Heaviside step function. In the definition of 
T the parameter is the resonant absorption coefficient defined such that 2a /7 is Beer's 
constant, where 7 is the halfwidth of the absorption line for the absorber. For this particular 
case we take 7 = A ph . The incident probability amplitude is assumed to be unity at the 
front face of the absorber (I = 0) for t — 0. According to Eq. (jJJ) the probability amplitude 
of the photon leaving a very thick sample (T > 1) in the forward direction is also close 
to unity for the leading pulse edge r ~ 0, the tail of the field being strongly absorbed and 
displaying an oscillatory behavior. This was experimentally detected for gamma photons in 
The process was called speed-up of the initial decay and dynamical beat 
for longer times (see, for example, the review Thus, the amplitude damping of the 



output radiation with increase of / (or T) does not follow Beer's law and the time integrated 



intensity of this field decreases as (2irT) x / 2 21 1 instead of exponential decrease e 2T . This 



is consistent with absorption in frequency domain of gamma quanta in a thick resonant 
absorber, which is calculated for a radiation source with a Lorentzian power spectrum J^. 

We argue below that the lack of absorption at r « results from the stepwise form of the 
initial part of the probability amplitude, i.e., it results from causality. In a cascade emission 
of two photons, the time at which the intermediate excited state e of the source particle 
is populated is known by detecting the first photon. From then on, and not earlier, the 
second photon can be emitted. Implicit in this scheme is the assumption that the emission 
of the first photon projects the atom into the intermediate level with absolute certainty. In 
quantum mechanical language, we project the two-photon state (produced in an atomic or 
nuclear cascade) onto a single photon state by detecting the first photon. This detection 
(the measurement) introduces a broad part to the otherwise narrow spectrum of the second 
photon characterized by the lifetime of the intermediate level r p h/2. A similar process in 
the photon gun of the first kind, i.e., the pumping of a single particle by a short laser pulse 
to the upper state followed by fast nonradiative decay to a fluorescent state, introduces a 
broad part to the spectrum of a single photon emitted in free space vacuum modes. This is 
the price paid for pulling the trigger at a desired time, i.e., for the knowledge of the time 
at which the excited state is populated. The time distribution of the detection probability 
of the photon emitted by the photon gun of the first kind looks similar to that observed by 
TDCM of photons emitted in an atomic or nuclear cascade (compare, for example, the plots 
for the photon probability presented in j^J and 14. lfi|). 



III. SPECTRALLY NARROW AND BROAD PARTS OF A SINGLE PHOTON 
EMITTED BY A SINGLE PARTICLE IN FREE SPACE VACUUM MODES 

In this section we give a brief outline of the methods used to treat the propagation of a 
single photon in a resonant absorber and introduce a decomposition of the photon spectrum 
(produced by a single particle in free space) in two components, one being a narrow and the 
other being broad. 

For simplicity we consider the probability amplitude of the emitted photon b(t) = 
b (t) exp(—iu t) normalized as b (t) = exp(— A ph t)0(t). The propagation of a single photon 
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in a thick absorber consisting of particles with the same resonant 



lifetime of the excited state as in the emitter is described by 



requency ujq and the same 



21] 



1 /-+0O 

bo(l,t) = — b (0 J v)e' iu< - t - l/c) - A ^ u)l dv, (5) 

27T 



where 



&o(0,z/) = - (6) 

A ph - iv 

is the Fourier transform of the incident photon bo(t) at the input I = and 

Ae 9 (v) = -^T- (7) 

7 — iv 

is the complex spectral function of the absorber. Here 7 is a half width of the absorption 
line of the individual particle and according to the imposed condition we have 7 = A p h- The 
parameter a = 27rN c n 2 u} /nc, defined in the previous section, contains the concentration 
N c of the resonant particles in the absorber, the refractive index n of the host in which 
these particles are incorporated, and the matrix element \x of the radiative transition e — g. 
The expression (J5J) was derived quantum mechanically by Harris , who considered three- 
dimensional resonant multiple scattering of a gamma photon in the forward direction by 
nuclei residing in a solid without regular structure of their positions. The same expression, 
Eq. (jSJ), also follows from the semiclassical approach considering the frequency dependence 
of the complex dielectric constant of the absorber 15]. Solving the Maxwell-Bloch equations 
or a small amplitude pulse with envelope bo(t) = exp(— A p h£)0(i) also gives the same result 



2l| . It was shown in Refs. [15J, |20|, |2l| that performing the integral in Eq. (JSJ) leads to Eq. 
©• 

To clarify the meaning of Eq. (j^J), we formally represent the Fourier transform of the 
single photon field, Eq. ©, as b (0, v) = b s (0, v) + fe a (0, v), where 



&a(M = X3T— (9) 



are symmetric and antisymmetric parts, respectively. The former has a Lorentzian shape 
and the latter resembles the dispersion part of the atomic response function. Their time 
domain counterparts are 

6 s (0,t) = ^exp(-A ph |t|), (10) 



|exp(-A ph t) if t>0, 
b a (0,t)=l if t = 0, (11) 

-±exp(A ph t) if t<0. 
The time dependence of b s (0, t) and b a (0,t) is shown in Fig. 2 by the dashed lines. We 
calculated the transmission of these components through a thick resonant absorber at the 
same conditions as for bo(T,t) with the help of Eq. (j3J), where the absorber length I is 
substituted by its effective thickness T = aol/j and 6 (0, ^) is substituted by b s (0, v) or 
b a (0, u). The result is 



b s (T,r) 



1 A ph r-T/2 
2 C 



if 
if 



r > 0, 
r < 0, 



(12) 



b a (T,r) 



§ (1 - e^/ 2 ) if 

_l e A ph r-T/2 if 



r > 0, 
r = 0, 
r < 0, 



(13) 



where 



^(T, r) = J (2 v^A^) ± J / e-^/ 2 J (2 y^A~^)dx, 



(14) 



and r = t — l/c is the local time at the output of the absorber of the physical length /. 
From these expressions one can find that for r = ±0 the probability amplitudes take values 
b s (T, ±0) = exp(— T/2)/2, b a (T, —0) = — exp(— T/2)/2, which are much smaller than 1 for 
large thickness (T ^> 1), and b a (T, +0) = 1 — exp(— T/2)/2, which is almost equal to the 
probability amplitude of the input radiation. For longer times (r > 0), the behavior of 
the probability amplitude of the antisymmetric part b a (T, r) is close to the total amplitude 
b (T,r), Eq. (J3J. The time evolution of the output probability amplitudes b s (T,r) and 
b a (T,r) for T = 10 is shown in Fig. 2 by the solid lines, (a) and (b) respectively. These 
plots clearly demonstrate that at the output of a thick sample the probability amplitude 
of the antisymmetric part becomes very close to the total probability amplitude of a single 
photon bo (T, r). We may conclude, that the narrow spectrum part of the photon is strongly 
absorbed in a thick sample while the broad component partially passes through and it 
defines the probability amplitude of the output radiation. The slight deviation of b a (T,r) 
from bo (T, r) is due to the reduced absorption of the far wings of the symmetric part of the 
photon spectrum in the case of A ph = 7. Its time domain counterpart b s (T,r) shows an 
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oscillatory deviation from the exponential behavior exp(— A p h |r| — T/2)/2 for r > 0. To be 
convinced that this deviation is due to the reduced absorption of the wings of the symmetric 
part of the photon spectrum, we consider the photon propagation through a medium with 
particles whose absorption line T is much broader than the spectral width of the input 
photon Aph- In this case we expect that the contribution of the symmetric part to the total 
output am plit ude bo(T,r) is negligible since it must be strongly absorbed. 

In Ref. |2l| it was shown that if, for example, we have inhomogeneous broadening of the 
absorption line in a macroscopic absorber, the output radiation field is described by Eq. (J5J) 
where the spectral function A eg {y) is substituted by 

AM = pSL. (16) 

with a total halfwidth V = 7 + 7^, where 7^ is inhomogeneous halfwidth. Generally, with 
the help of the linear response approximation, one can derive the same expression for a weak 
radiation field transmitted through a resonant absorber having an absorption halfwidth T 
(natural, homogeneously or inhomogeneously broadened by the neighbors of the resonant 
particles in the absorber), which is different from the spectral halfwidth of the incoming 
radiation A p h (A p h <C T). We calculated the transmission of the probability amplitudes 
of the symmetric and antisymmetric parts of incoming radiation for such a sample with 
spectral function Ar(u). The result is 

I e -A ph r-T_ + j n ) if > o 
Ul,r)={ 2 i [ ( ' - (16) 

l e -A ph .-T_ +/+(/)T) if r>Q; 

b a (l,r)={ |(l-e" T+ ) if r = 0, (17) 

_I e A ph r-T + if T < 0> 

f±(l,T) = ^[g(T_,T)±g(T + ,T)], (18) 

9(T±, t) = J ± e- (T± - x) J (2^x(T ± A ph )r)dx. (19) 

Here T± = oiqI/(T ± A p h) is the effective thickness of the absorber. Formally, with the 
decrease of A ph (A ph — > 0) the function /_(/,r) in the solution b s (l,r) tends to zero since 



where 



its components g(T + ,r) and g(T-,r) nearly cancel each other. However, their difference is 
larger than the first term (1/2) exp(— A p hT — TL) in b s (l, r) for r > 0, and, hence, we have 
residual oscillations of the symmetric part of the photon amplitude for r > 0. In spite of 
this, for r > 0, it is the antisymmetric part b a (l,r) that gives the dominant contribution 
to the total probability amplitude of the output photon. This amplitude, b (l,r), can be 
approximated by the two main terms obtained from the successive integration by parts of 
the integrals in Eqs. (fTBj) . (fTTj) . which are 



b {l,r) = e 



-Ft 



J {2^a lT) + (r - A ph )r — 



9(t). (20) 



J\(x) is the Bessel function of the first order. It is remarkable that the first and also main 
term of this expression does not depend on the spectral width A p h of the incoming photon. 
Thus, for Aph <C T the amplitude of the output radiation is mainly defined by the parameters 
of the sample and not those characterizing the input photon. 

The symmetric part of the input photon spectrum has a Lorentzian shape with halfwidth 
Aph that is narrow with respect to the absorption linewidth, i.e., A p h <C T. Therefore, this 
part is strongly absorbed in a thick sample. The antisymmetric part partially passes through 
a thick sample and what comes out is due to the far wings of b a (0, u), which are proportional 
to 1/u and independent on A ph . This is the reason why the output radiation does not contain 
information about A p h. Therefore, formally we refer to this part of the photon spectrum 
as broad. Fig. 3(a) shows a comparison of the time dependence of the antisymmetric part 
b a (l) T ) with the numerically calculated integral (0) for the total amplitude b (l, r), where the 
spectral function is Ay{v). These plots show that the probability amplitude of the output 
radiation essentially coincides with the probability amplitude of the antisymmetric part 
b a (l,r), confirming the heuristic argument that the symmetric part is strongly absorbed. 

Another argument supporting the concept of the broad spectrum component of the pho- 
ton can be given with the help of the time integrated intensity of the output radiation. For a 
classical field, the total energy transmitted through a unit area at distance I is proportional 
to 

/+oo 
b(l,r)b*(l,T)dr. (21) 
-oo 

For a single photon, this value is actually proportional to the number of counts of the second 
detector in a wide time window obtained without TDCM, i.e., without the first detector in 
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the scheme Fig. 1(b). Eq. (|21|) can be transformed as 

1 p+oo 

U ( l ) = 2^J Kl,v)b*(l,v)du, (22) 

where b(l, v) = b(0, v) exp[ — ^4 m (z/)/] and A m (z/) is A. e g{v) or Ap(z/). 

In the case of a macroscopic absorber with a narrow spectral function A eg {y) of halfwidth 
A p h, we have for the symmetric, b s (l,v), and antisymmetric, b a {l,v), parts of the output 
photon: 

U S {T) = e- T [I (T) - h(T)]U (0)/2, (23) 



U a {T) = e- T [I (T) + J 1 (T)]f/ (0)/2, (24) 

where U (0) = r ph /2 is the time integrated intensity of the input field b (0, r), T = a //Ap h , 
Io{T) and I\{T) are the modified Bessel functions of zero and first order, respectively. For the 
time integrated intensity of the total field b (l, v) = b s (l, v) + b a (l, u), we have the expression 
JUT) = U S {T) + U a {T), which is consistent with U (T) = exp(-T)I {T)U (0) found in Ref. 
2l|. 

For an absorber with a broad spectral function Ar(z/) with a halfwidth T > A p h, we have 



U a (T b ) = U + {T b ) 



A 



P h 



r 



[Ui{T b ) — U 2 (T b )} , 



(25) 



where 



U a (T b ) = U-(T b 



A 



P h 



U±(T b 



-2aT h 



1 ±4a 2 



A 



ph 



r 



T b 



U 2 (T b ), 



U (0) 



(26) 



(27) 



U x (T b ) = 2a 2 U (0) / e- 2a{Tb - x) - x I (x)dx, 
Jo 



U 2 (T b ) = 4a 3 f/ (0) [ T \T b - x)e- 2 < T ^~ x T (x)di 
Jo 

a = T 2 / (r 2 — A 2 h ) and T b = a l/T is the effective thickness of the absorber at the center of 
the broad absorption line. Fig. Efb) shows the thickness dependence of the time integrated 
intensity of the symmetric, U s (T b ), (dash-dot line) and antisymmetric, U a (T b ), (solid line) 
parts of a single photon. The antisymmetric part does not follow Beer's law, while the 



(28) 



(29) 



11 



symmetric part follows this law up to a certain thickness (T& < 5) and its contribution to 
Uq(T) becomes negligible for T > 1.5. For an absorber with a narrow absorption line of a 
halfwidth A p h, both parts (not shown on the plots) do not follow Beer's law for any distance 
and their T-dependence resembles that for U a (Tb) shown by thick solid line in Fig. 3(b). 

If A m {y) = A eg {y), Eq. (J22j) for U(l) coincides with a formula for the transmission of 
gamma-photons through a thick absorber of length I (exact resonance case if the central fre- 
quency of the photon equals the resonant frequency of an absorber), well known in frequency 
domain Mossbauer spectroscopy 23j. This means that the Lorentzian power spectrum of 
gamma photons, adopted in the general formula j^j], implies causality, since the relevant 
spectrum of the photon probability amplitude must satisfy Eq. (JTJ) or Eq. ©, both assum- 
ing a time-asymmetric photon envelope. If the photon spectrum would be defined only by 
the symmetric part 6 s (0,z/), Eq. (jHJ), the power spectrum of such a photon is Lorentzian 
squared and all peculiarities mentioned above (violation of Beer's law etc.) would not be 
present. 

Summarizing we conclude that the spectrum of the photon emitted in free space consists 
of two parts: a narrow and a broad part. The narrow part comes from the symmetric 
component b s (0,t) of the photon probability amplitude. The amplitude of this part decays 
with propagation distance and its absorption follows Beer's law in a medium with a wide 
absorption line. In a medium with a narrow absorption line, comparable with the spectral 
width of the symmetric part, this component is also strongly absorbed. But due to the 
reduced absorption of its far spectrum wings, the frequency integrated absorption of the 
symmetric component does not follow Beer's law. The antisymmetric part of the photon 
always violates Beer's law, its amplitude at r = +0 is nearly 1 for any thick absorber, then 
this part decays with the rate defined by the spectral width of the absorber and by its 
effective thickness. 

In the next Section we propose to filter out the broad part of the photon spectrum with 
the help of electromagnetically induced transparency. The same mechanism can also reduce 
the intrinsic spectrum width of the source photon, i.e., the width of the narrow part. 
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IV. PHOTON FILTERING THROUGH THE EIT WINDOW 



Electromagnetically induced transparency (EIT) allows transmission of the radiation field 
through a thick sample without absorption and with appreciable delay due to the reduced 



group velocity Ref. 



24j . Recently it was shown that the EIT window allows also a separation 



of the narrow and broad components of the radiation field in time |25| . Here we apply these 
results to separate the narrow and broad components of a single photon emitted in free 
space. For simplicity we consider a resonant absorber described by a three-level scheme 
with a DC coupling field shown in Fig. 4. This scheme was experimentally realized for 
atomic hydrogen whose excited electron states 2p and 2s (degenerate in energy) are coupled 
by an electric DC field I26]. State 2s is a metastable state with a long lifetime. Therefore, 
the coupling creates a narrow transparency window at the center of the absorption line. 

For a weak input radiation field being in exact resonance with e — g transition, the output 
field from the absorber with EIT window is described by Eq. (jSJ) with the spectral function 
Agg^) replaced by |2q] 

= (30) 

(r - iu)^y m ~ w) + YL 2 

where T is the halfwidth of the unperturbed absorption line for the transition from the 
ground state g to the excited state e, Q is the coupling parameter defined by the strength 
of the DC field, which couples the excited state e with the metastable state m (see Fig. 4), 
and 7 m is the half decay rate of the metastable state m. The EIT hole is deep if the coupling 
parameter satisfies the condition Q 2 ^> 7 m r and its halfwidth is approximately A e i t = Q 2 /F. 

Following Ref. 2^1, we apply the adiabatic expansion of the spectral function A c \ t (v) 
near the center of the EIT hole, v — 0, for the approximate calculation of the integral in Eq. 
(JSJ) describing the output photon probability amplitude. This method is valid if the EIT 
hole is present at the line center, i.e., if the condition Q 2 > 7 m r is satisfied. The method 
gives a nice approximation of the spectrally narrow, adiabatic part of the output radiation 
filtered through the EIT window in a thick absorber. The nonadiabatic, spectrally broad 
part of the output radiation will be treated separately, as in (2^. In the adiabatic expansion 
it is sufficient to retain the three main terms (see Section III in Ref. 

A eit (is)l^T eit -ivt d + v 2 /A 2 cS , (31) 
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containing three important parameters of the EIT medium, i.e., 



T <* - r wb (32) 



T ' r (TOF (33) 



*cff 



+ ImTf 



r 6 r[fi2(r + 2 Tm )-73] 

The first parameter, T e j t , is the EIT reduced effective thickness of the absorber, which is 
defined by the residual absorption at the bottom of the EIT hole [y = 0). Here, we define the 
effective thickness with respect to the absorption without the coupling Q, i.e., T& = aol/T. 
We have T eit /T b pa /S. eit <C 1 if fi 2 ^> 7 m r. The second parameter, approximated as 
id ~ 7l/Aeit, is the delay time of the radiation field associated with the slow group velocity 
V g = c/(l + ctd/l). It is defined by the slope of the steep dispersion at the center of the 
EIT window. The third parameter, A e fj, is the effective width of the EIT window for a 
thick sample, which narrows with distance as ~ 1/y/Ti. This parameter is approximated 
as A e fj pa A^t/y/Tl and it is responsible for the pulse broadening in time or its spectrum 
narrowing with distance. With A e i t (v), given by Eq. (j31j) . the integral in Eq.fjSJ), which can 
be considered as the inverse Fourier transform, is calculated with the help of the convolution 
theorem. This gives an integral representation of the adiabatic solution for the probability 
amplitude of the output photon 

b 0A (l, t) = =±I b 01 (l, n - t d )e^ A ^-^ 2 dn. (35) 

The exponential part of the integrand takes into account the third term of the adiabatic ex- 
pansion (Eq. (JHU). The other part of the integrand, b i(l, T\—t d ) = b (0, T\ — td) exp(— T eit ), 
is the probability amplitude obtained using only the first two terms of the adiabatic expan- 
sion, i.e., Ae it (v)l pa T eit — ivt d . Explicitly this function is 

boi(l,r-t d ) =e- A ^-^- T ^Q(r-t d ). (36) 

The integral in Eq. (|33j) can be calculated and gives 

b 0A {l, r) = 0+(Z, T )e- T "- A *<r- t *\ (37) 
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where 



Ml, t) = \l ± erf ( ^(r - t d ) T 




A ph 



Aeff 



) 



(38) 



erf(x) is the Error function, which is defined as an odd function, i.e., erf(— x) = — erf(x). 
The function 0_(/,r) will be used to describe the probability amplitude of the symmetric 
and antisymmetric parts of the photon. 

If A ph <C A cff , the function </>+(/, r) is almost unity for r t d + 2A ( J ff 1 , and it is almost 
zero if r <C td — 2A ( J ff 1 . In the transient domain, i d — 2A cff 1 < r < t d + 2A^ ff 1 , the function 
<p + (l,r) gradually increases from to 1. Fig. 5 shows the time dependence of the function 
r) for Aph/A e fj = 0.1 (solid line) and for A p h/A e ff = (dotted line). We see that if the 
spectral width of the narrow part of the photon, 2A p h, is much smaller than the effective 
width, 2A fj, of the EIT window for a thick sample (T& ^> 1), the exponential tail of the 
delayed photon is almost not affected by the EIT hole, while the front edge is smoothened 
and spread around td in the domain td — 2A cff 1 < r < t d + 2A^ ff 1 . Thus, the narrow part 
is not affected in this case while the broad part is filtered out. Actually the broad part is 
present in the output signal, but it is not delayed, as will be shown below. 

If the spectral width of the narrow part of the photon is much larger than the transparency 
window, A ph ^> A eff , then the function 6 01 (/,r 1 — td) in Eq. (|33j) can be approximated by 
a delta function, 5{r\ — td), and its convolution with the Gaussian function is 6oa(^ t ) ~ 
exp [— iA^ ff (r — td) 2 ]- Thus, for A ph ^> A eff the photon envelope will spread in time and 
narrow in spectrum acquiring a Gaussian shape with a temporal width 4/A e ff. 

Such time-broadening of the pulse almost conserves its area 



whose value is affected only by the residual absorption defined by T eit , which can be very 
small (T eit <C 1)- Here 9a(0) = 2/i/A p h is the area of the input pulse b (0,t). The pulse 
area is a quantitative parameter specifying the field interaction with a two-level atom. For 
example, the probability amplitudes of ground and excited states of a two-level atom inter- 
acting with a short resonant pulse (shorter than any relaxation time) are given by cos(6 l /2) 
and ism(9/2) respectively, where 9 is a pulse area |27l |. 

The time integrated value of the pulse intensity ~ |&oa(^ t )| 2 , which is essentially the 



9 A (l) = / 2fib 0A (l, r)dt = 9 A (0)e- T ^\ 



(39) 



J — oo 
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average detection probability of the photon, 



U A (l)= / \b A(l,r)\ 2 dr 



(40) 



reduces with distance due to broadening of its envelope as 



U A (l) = U (0)e 



2T eit +2A2 




(41) 



If A p h -C A e ff/ a/2, the output energy does not change appreciably due to the pulse broad- 
ening in time. If A p h ^> A c ff/ y/2, the energy is reduced by a factor that is the ratio of the 
EIT hole width and the spectral width of the photon ~ A e fj/A p h. 

So far, we have discussed the propagation of the adiabatic component of the probabil- 
ity amplitude of the photon whose spectrum is confined in the EIT window. This is the 
spectrally narrow part. The spectrally broad part is not sensitive to the narrow frequency 
domain of the absorption spectrum around v = 0, where the EIT hole is located. The time 
evolution of this spectrally broad part of the photon in a thick resonant absorber is mostly 
defined by the far wings of the absorption line, which decrease as 1/ z/ 2 , and by the dispersion 
component, which drops even more slowly ~ 1/u. Therefore, to describe the propagation of 
the spectrally broad part we neglect in the function A(y) the coupling Q, which produces 
the EIT hole and does not affect the wings of the absorption line. Then, the broad, nonadi- 
abatic component of the photon probability amplitude, bo(l,r), is governed by the integral 
© with the spectral function A^{y) given by Eq. (fI5jl. 

We consider two cases. If /S. p h = T, we have the case of equal spectral widths of the 
incoming photon and of the absorption line for a single particle. Then the nonadiabatic 
broad part is &at(Z,t) = boN(l, T ) exp(— zcg> t), where &ojv(/,t) is described by Eq. (J2J). The 
second case corresponds to the situation where the photon spectral width is smaller than the 
absorption linewidth (A ph < T). Then &oat(£, t ) is given by the sum of b s (l, r), Eq. (|TB|). and 
b a (l,r), Eq. (fT7|). or by the simplified expression (|2Uj) . In both cases the output probability 
amplitude of the photon is the sum of the adiabatic (spectrally narrow) and non-adiabatic 
(spectrally broad) amplitudes 



First we consider the case A ph = j m . Then the two terms, A ph t rf and T eit , are almost 
equal and they nearly cancel each other in Eq. (p?7J). reducing the adiabatic part of the 



btat(l, r) = b 0A (l, t) + b 0N (l, t). 



(42) 
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photon amplitude 6 a(^,t) to 

M/,r)^<MZ,r) e - A p^. (43) 

Figure 6(a) shows the time dependence of the photon probability amplitude b tot (l,r) (bold 
line). This is our analytical solution, Eq. (J4*2|) . It coincides almost perfectly with the 
numerical evaluation of Eq. (jSJ) with y4 eit (z/), which is an exact expression. The effective 
thickness is taken as T = 30, the other parameters are T = 107 m and Q = 2T. The 
probability amplitude of the photon at the same distance with no resonant absorber between 
the source and the detector, 6 (0,t), is shown for comparison (thin line). The dashed line 
shows the probability amplitude boi(l,r — td), Eq. (}3l)|) . which does not take into account 
the spectrum narrowing of the photon. The pulse area of boi(l,r — td) (the area under 
the dashed line) coincides with the area of the adiabatic part, b A(l,T). Observing only 
the tail of the photon envelope, one cannot distinguish the delayed photon from the one 
that would reach the detector without absorber. However, the leading edge of the photon 
envelope is split in two parts. The non-adiabatic part is fast with a duration proportional 
to 1/ (T&r) = l/(a /) and its amplitude equals unity at r = 0. This part originates from the 
spectrally broad component of the photon. The thicker the sample, the shorter this part 
will be. The adiabatic component is slow and propagates with group velocity V g . Its leading 
edge is broadened in time from t = td — 2A^ to r = td + SA^ 1 . The slow part of the photon 
does not change appreciably in shape except the smoothening of the leading edge. This is 
because for the chosen values of the parameters T, Q, and T the effective width of the EIT 
hole is larger than the width of the photon spectrum, i.e., A e g > A p h- For our numerical 
example A e fj/A p h = 6.9. The delay time of the photon, td, is 1.4 times longer than the 
lifetime of the excited state of the source particle Tiif e = r p h/2. 

Since A c fj « Q 2 /T\^Tl, one can decrease the effective width of the hole and make it 
narrower than the spectral width of the photon, A e g < A p h, by reducing the coupling Q or 
increasing the thickness I of the sample and hence increasing the effective thickness T&. Then, 
it can be shown that the shape of the slow photon amplitude reduces to almost a Gaussian. 
However, its far tail for long times cannot decay more slowly than exp(— j m t). Thus, if 
A p h = 7 m and an EIT hole is present, the far tail of the photon probability amplitude 
is always indistinguishable from the case without absorber. This is because the induced 
coherence of the states g and m has the same lifetime as the lifetime of the photon probability 
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amplitude, both being caused, for example, by the spontaneous decay of the excited (for the 
source) and metastable (for the absorber) states. The g—m coherence (i.e., an antisymmetric 
superposition of states g and m) is at the origin of destructive interference. It can be shown 
that this coherence reduces the absorption and as such supports the photon propagation 
without absorption. 

The indistinguishability of the tails of the slow photon and a photon that does not 
propagate through an absorber can be removed if A p h ^> 7. In this case the slow part of 
the photon envelope transforms to an almost Gaussian shape. Most part of its amplitude 
exceeds the probability amplitude of the photon propagating without absorber. The area 
under its envelope coincides with the area of b i(l,r — to), i.e., with 6a(1) given by Eq. 
()39|) . Fig. 6(b) shows the time dependence of the photon probability amplitude 6 tot (/,r) if 
Aph = T (bold line). The other parameters are the same as for Fig. 6(a). Comparing the 
fast, nonadiabatic parts of the probability amplitudes in Fig. 6(a) and 6(b), we find that 
they are almost the same in spite of the difference of the halfwidth of the photon spectrum 
for (a) and (b), being 7 m and T, respectively. This is because the spectrally broad component 
of the photon, responsible for the nonadiabatic contribution, originates from the stepwise 
rise of the photon probability amplitude at t — and does not depend on the spectral width 
of the radiation source. 

Concluding this section we analyze the evolution of the symmetric and antisymmetric 
components of the photon spectrum with distance in the EIT medium. Their adiabatic 
6 s (A)(^,t), 6 a (A)(/,r) and nonadiabatic b s (N)(l,T~), b a ^)(l,r) counterparts can be calculated 
similarly to b A(l, r) and box (I, r). Applying the same procedure as before we obtain 



where the adiabatic symmetric b s (A)(l,T) and antisymmetric b a (A){l,T) counterparts are 



(44) 



b s{A) (l,T) = [R + (l,T)+R„(l,T)}/2, 



(45) 



b a(A) (l,T) = [R + (l,T)-R_(l,T)]/2. 



(46) 



The function R±(l, r) is defined as 



R±(l, r) = (f)±(l, t) exp [-T eit =f A ph (r - t d )) . 



(47) 
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The nonadiabatic parts & s (jv)(/, t), b a ^(l,r) are calculated in Sec. III. They are defined by 
the output field for the symmetric b s (l,r) and antisymmetric b a (l,r) components given by 
Eqs. (|12j ) .()13 )) for A ph = T and by Eqs. (fTB |) .(|T7 |) for A ph = j m . For the latter case, 7 is 
substituted by / ~f m in Eqs. (|16j) . (jl7|) . Fig. 7 shows the time dependence of the symmetric 
(a) and antisymmetric (b) components of the photon probability amplitudes with a narrow 
spectrum, A p h = j m . These plots clearly demonstrate that the broad component of the 
photon spectrum is present only in the antisymmetric time domain counterpart. Both 
counterparts are delayed due to the reduced group velocity and they are smoothened. The 
discontinuities in the time dependence of the time derivative for the symmetric part and of 
the amplitude for the antisymmetric part are removed. 



V. APPLICATION OF THE EIT FILTERING 



The construction of a quantum network, consisting of quantum nodes and interconnecting 
channels, is a challenge in quantum communication and computation science. It is natural 
to use as nodes matter in the form of individual atoms or atomic ensembles. Storage of a 



quantum state of light in an EIT medium 



28j or in an extended ensemble of atoms with 



the photon-echo technique J29} looks very attractive. The experimental realization of such 
protocols with a classical radiation field |30| and faint laser pulses imitating a single photon 
radiation field |31] showed the feasibility of storing and retrieval of a radiation field and 
demonstrated a delayed self-interference for a single photon wave packet in matter. If a 
single photon source of the first kind is used in these protocols, the fidelity of the mapping 
of a quantum state of light to matter can be quite low. To show this we introduce two 
quantum states 

1 j, \ '/«r A ph exp(-ik- r ) , . . . 

\b a ) = -zV2 2^ffk p 2 A2 |lk>, (48) 



and 



1 1 \ /ttv^ z/ fc exp(-«k- r ) , , , 
\b a ) = V22^g k |l k ), (49) 

k u k -r ^ P h 

which correspond to the symmetric and antisymmetric parts of the photon emitted in free 
space. The states are orthogonal and normalized: (b s \b s ) = (b a \b a ) = 1; (b a \b s ) = 0. This 
can be easily verified taking into account that the Fock states |lk) are orthogonal and the 
sum over k is replaced by an integral over all frequency modes. The orthogonality of the 
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states comes from the fact that the symmetric part of a photon, \b s ), is an even function of 
Vk while the antisymmetric part, |6 a ), is an odd function of v^. 

Formally we can represent the radiation state \b), emitted by a single quantum particle 
in free space, [see Eq. (JTJ) in Sec. II] as 

\b) = ±= \b s ) + ^ |6a> • (50) 

Such an expression allows us to make a strong statement about a single photon interacting 
with a single atom or an ensemble of atoms. First, it is evident that the photon emitted by 
a single particle populates equally both states, \b s ) and \b a ). In Sec. Ill, we showed that the 
symmetric, \b s ), and antisymmetric, \b s ), parts of the radiation field interact differently with 
resonant atoms (quantum absorbers in general). The antisymmetric part almost does not 
interact with an atom while the symmetric part does. Also, it is known (see, for example, 
Ref. H) 

that the probability of excitation of an atom by a radiation field with a temporal 
duration much shorter than the decay time of the atomic coherence is defined by the pulse 
area. This area is zero for the antisymmetric part of the photon and it has a finite value 
for the symmetric part. Thus, a photon emitted by a photon gun of the first kind has 
only a 50% chance to interact with a target atom and 50% chance to miss the target atom. 
Such a poor score would make the use of a single photon source of the first kind in quantum 
computing or information storage unreliable. An experimentalist performing a measurement 
cannot distinguish whether missing a target is the result of quantum statistics (i.e., the result 
of a small interaction probability) or a purely quantum mechanical result originating from 
causality (i.e., the result of pulling the trigger, which leads to populating the noninteracting 
state \b a )). Therefore, a qubit preparation in the form of matter excitation or any other 
operation with such a qubit using a single photon source of the first kind would become very 
uncertain and hence unreliable. 

However, filtering a single photon through an EIT window may help to make it much 
more reliable for quantum computing and information storage. At the output of a thick EIT 
filter, the broadband, \b a ), and narrowband, \b s ), parts of the photon are well separated in 
time. The broadband part can be removed by time gating or a shutter and, for example, sent 
to an auxiliary channel. Then, the probability amplitude of the photon at the output of the 
EIT filter becomes bell-shaped (see Sec. IV). Sending the removed broadband part of the 
photon to a detector allows to purify the photon state. If the detector does not " click" , all 
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probability amplitude is collected in a state that is a wave packet with a Gaussian envelope 



If it clicks, no photon with a Gaussian shape is present. Such a detector in a purification 
scheme can be omitted, since what comes out of the EIT filter with appropriate time gating 
is always the pure state (|51|). However, this auxiliary detector may help to conclude that 
the photon, emitted on demand, failed to pass through the EIT filter and we have to repeat 
the operation. 

One could argue that filtering the radiation emitted by a single photon source reduces the 
total probability amplitude of the output photon while the interaction probability of such a 
photon with a target atom (or atoms) remains almost the same. So, what would be the gain 
and what would be the advantage of the EIT filtering? The gain comes from the removal 
of the broadband part, which does not interact with the target atom and hence produces 
a count at the detector placed behind the target. As such it is a false count carrying no 
information since the photon is assumed to be "stored" in the atom (atoms), but it escapes 
the atom-field interaction due to the |6 a )-component. Therefore, the false counts can be 
considered as noise in an information storage process. When applying filtering, this noise 
would reduce to zero, while with no filtering we have to compare the atom-field interaction 
probability with the quite high detection probability (~50%) given by the broad spectral 
component of the photon \b a ). 

A further narrowing of the symmetric part of the photon in the case A p h > A eit works in 
the same way. We lose the "brightness" of the source or the detection probability, but gain 
in spectral resolution. Actually we do not lose the "useful" energy if our aim is to improve 
the spectral resolution by sharpening the line associated to the source photon, i.e., to make 
it more selective in the excitation of a particular target atom. Selective excitation takes 
place if the spectral width of the source photon equals the spectral width of the absorption 
line of the selected atom and if their frequencies are in resonance. If our filter is designed 
such that these requirements for the output radiation are fulfilled, only the "useful" spectral 
content of the source photon is transmitted and the rest is suppressed. As a result, the 
target atom interacts with such a photon with the same probability as in the case of no filter 
and the "useful" brightness of the source does not change. Otherwise, with no filter the 
inherent broad part of the photon spectrum not interacting with the target atom (atoms) 




(51) 
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would produce a click at the detector placed behind the target showing that the atom-field 
interaction failed to happen. 

To conclude this section, we consider the interaction of a Gaussian shape photon with an 
ensemble of resonant atoms. We define the Gaussian wave packet as 



IM = ET^^ k - r °^ /A ^|l k ), (52) 

normalized such that it contains only one photon. The associated single photon field 
&g(*) — (0| E ( + )(r, t) \ba) is &g(0, t) = b 0G (0, t) exp(— iu t), where for simplicity its am- 
plitude &og(0, t) = exp(— A 2 h t 2 /4) is normalized to unity for t = 0. The propagation of such 
a photon in a thick absorber consisting of particles with the same resonant frequency uj 
is described by Eq. (jSJ). We consider the case as in Sec. Ill, if the absorption line in the 
absorber is inhomogeneously broadened with half width T [see Eq. (fT3j) ]. To compare the 
absorption of a single photon emitted by a source of the first kind with that for a photon 
with a Gaussian envelope, we address the case of a narrow spectral width of a single photon: 
Apj, <C T. In this case the integral Eq.fjHJ) for a photon with the Gaussian envelope, which 
describes its transmission through an absorptive medium, can be calculated expanding the 
function v4r(z/), Eq. (fT5j) . in a power series of v near v = 0. Keeping only three terms of the 
expansion, we have 

£ (!+£-£). (53) 
Then the integral is calculated analytically, i.e., 



b (T,r) « rjexp 



-T - V 2 A 2 ph ( r 



rj-\ \ 2 



(54) 



where T = aol/T, r = t — l/c is the local time, r\ = l/y/l — JT, and / = (A p h/T) 2 . This 
approximation is valid if f'T < 1. The transmission function for a Gaussian photon, Ug{T), 
can be easily calculated 

U(T) = ^e~ 2T . (55) 

For r] ~ 1, the deviation from the Beer's law is negligible, which shows a strong atom-field 
interaction for such a photon. Thus, the state of a single photon with a Gaussian envelope 
can be mapped into an atomic ensemble. Recently, the mapping and retrieval of a classical 
: ield state in an ensemble of resonant impurities in a solid was experimentally demonstrated 

ai. 
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VI. LEVEL MIXING INDUCED TRANSPARENCY FOR GAMMA RADIATION 



A similar analysis is applicable in t he g amma domain, in particular to reexamine the 
recent level crossing experiments in 57 Fe 13, [mJ]. In these references by means of Mossbauer 
spectroscopy [3| the absorption of the 14.4 keV photons, emitted by a radioactive single line 
source 57 Co (CoRh), has been experimentally studied in the natural mineral siderite (FeCOs) 
containing 57 Fe nuclei. At the crossing of two hyperfme levels mi = —3/2 and m 2 = +1/2 
in the first excited state of 57 Fe j^], an appreciable deficit of gamma-photon absorption 
was found. The experiments were performed in frequency domain by Doppler shifting the 
frequency of the radiation source and detecting the 14.4 keV photon, which passed through 
a thick resonant absorber containing 57 Fe. One can expect that time domain experiments 
with TDCM could show a photon delay at the level crossing of 57 Fe. We assume that the 
fluctuating electron spin of Fe 2+ in FeC03 is the dominant source of the line broadening 
for the transition g — e in which the ground state level g is m g = —1/2 and the excited 
state level e is mi = —3/2. Our assumption is also supported by studies in Ref. jsfj. 
Meanwhile, the transition g — to, where m is the excited state level m 2 = +1/2 has mostly a 
natural broadening caused by spontaneous decay via the emission of a gamma photon (into 
47r angle) or via electron conversion [3^. If the excited levels mi = —3/2 and m 2 = +1/2 
are coupled by symmetry breaking interaction, then we have an EIT scheme with a coupling 
of two excited levels having different line broadening mechanisms, one is due to the electron 
spin fluctuations producing the linewidth 2r and the other due to natural broadening, 27™, 
so that T > 7 m . This EIT scheme will produce a slow photon shown in Fig. 6(a), where 
A ph = 7m . Aether polity for EIT is to apply an rf mixing of two .evels in the excited 
state of 57 Fe |37|. 

The group velocity of the single photon wave-packet can be reduced many times inside an 
absorber due to the steep dispersion in the EIT window. If the propagation time of a slow 
photon in the absorber is much longer than the mean dwell time between successive photons 
coming from the source, one can collect several photons in the sample and then release them 
when required by means of a sudden removal of all nuclei from resonance, producing a short 
burst of radiation. 
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VII. DISCUSSION 



We have shown that a photon emitted by a single particle in free space contains spectrally 
narrow and broad components with equal probability amplitudes. The broad component 
has very low absorption in an ensemble of two-level absorbers or a single atom. Therefore, 
a mapping of the state of such a photon in the form of matter excitation has ~50% proba- 
bility, which is low for any realistic protocol of quantum information storage and quantum 
computing. This is very different from what can be expected for a single photon with a 
Gaussian time-envelope, which demonstrates a high probability of atom-field interaction. 

Filtering out the broad component of the photon with the help of an EIT medium may 
produce a single photon that is (i) spectrally narrow and (ii) of Gaussian shape (if A ph > 
A e it) in spite of an initially highly asymmetric time dependence of the probability amplitude 
of the photon coming from the source. Photon reshaping can be made at low cost in time 
and energy. The narrow and broad spectral components of the photon are separated in time 
and space because of different group velocities. The narrow part is broadened in time (or 
spectrally narrowed) and delayed, while the broad part transforms into a fast photon with a 
very short duration. The separation of these components, for example by time gating, may 
produce a photon with a much narrower spectrum, even with respect to the narrow part A p h 
if Aph > A eit . The spectrum narrowing of the narrow part takes place without appreciable 
amplitude probability loss. This is because the time integrated probability amplitude of the 
photon, which is analogous to the classical pulse area, is conserved and is not affected by the 
spectrum narrowing of the photon. This area is just a parameter quantifying the atom-field 
interaction. However, the time integrated probability of the photon (its amplitude squared), 
which is analogous to the energy of a classical pulse, decreases inversely proportional to the 
square root of the effective thickness of the absorber. This decrease is much smaller than 
the value given by Beer's law for monochromatic radiation. 

Recently, another type of single photon sources based on spontaneous Raman scattering of 
a laser field in an extended medium was reported j^-j^]. They are designed to implement 
the DLCZ protocol |43( for long-distance quantum communication with atomic ensembles 
and linear optics . Spontaneous Raman scattering of a low intensity laser field in an extended 
ensemble of three-level atoms with states \g), \s), |e), is capable to produce a low intensity 
Stokes field in the forward direction with an average photon number smaller than unity. In 
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this protocol, the laser field is strongly detuned from the resonant transition \g) — > |e) to 
insure low population of the excited state atoms, |e). If the atoms are initially prepared in 
the ground state \g), which the laser field excites, the spin-wave (coherence of the ground 
states | g) and |s) distributed wave-like in the ensemble of three-level atoms) is excited along 
with the scattering of one Raman photon into a wave packet propagating in the forward 
direction. This event is probabilistic and its success heralds the preparation of one quantum 
spin-wave excitation stored in the atomic ensemble. It also demands the same detuning 
of the Stokes photon from resonance as for the laser field such that the Raman resonance 
condition for the two-quantum process \g) — > |e) — > \s) is satisfied, i.e., the frequency 
difference of the laser field and the scattered radiation equals the frequency of the transition 
between the ground states \g) — > \s). If the Raman resonance condition is not satisfied, 
with high probability no spin- wave but rather the population of the ground state \s) for a 
single atom (randomly localized spin excitation) in the atomic ensemble is produced after 
spontaneous Raman scattering on this atom with emission of a photon with a resonant 
frequency u es = [E e — E s )/h, where E e and E s are the energies of states |e) and \s). 
Therefore, this kind of a single photon source satisfying the Raman resonance condition 
needs a post-selection procedure. It implies the generation of pairs of single photons, Stokes 
and anti-Stokes. One photon of the pair (Stokes) is generated by the Raman scattering 
described above. The laser producing this photon is named a "write laser". It generates 
the photon and one quantum of the spin wave excitation. Another photon of the pair (anti- 
Stokes) is generated by a "read laser". Depending on the chosen scheme, it is applied to the 
transition \s) — > |e) (in a three-level scheme) or \s) — > |e') (in a four-level scheme), where 
|e') is an other excited state of the atom. The read laser generates an anti-Stokes photon 
in a two-quantum process \s) —>■ \e) —>■ \g) or \s) — > \e') — > \g). Detection of a single anti- 
Stokes photon after the read laser pulse (post-selection) warrants the presence of a single 
quantum of a spin-wave created by the write pulse and the generation of a single photon 
state in the Stokes field satisfying the Raman resonance condition. To reduce two-photon 
scattering in the Stokes mode, the probability of a single photon emission, n, must be small 
(n < 1). Then the two-photon emission probability of the Stokes field becomes even smaller 
~ n 2 . In this way a correlated pair of single photon Stokes and anti-Stokes fields is created. 
Stokes field generation in an atomic ensemble has collectively enhanced the coupling to a 
certain optical mode due to many-atom interference effects [44]. This mode is defined by 
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Raman photon scattering into a wave packet propagating forward along the write laser field 
and it is supported by the spin-wave if the Raman resonance condition is satisfied. All 
Raman scattering trajectories of the Stokes photon interfere constructively in the forward 
direction. This is very similar to the enhancement of single photon scattering in the forward 
direction in an ensemble of two-level atoms considered in Sec. II of our paper. The difference 
comes from the spectral content of the produced photon. Due to the large detuning of the 
write laser pulse and the post-selection of the Stokes photon satisfying the Raman resonance 
condition, which warrants a single quantum spin-wave excitation, the time envelope of such 
a photon is completely defined by the spectral properties of the write pulse and not by the 
lifetime of the excited state |e). This was clearly shown for the Stokes photon with photon 
number n < 1 by Lukin's group in Ref. |4(3] . For n > 1, the spontaneous generation of the 
Stokes field changes to stimulated emission, introducing a time-asymmetry of the scattered 
field envelope. To reduce the spontaneous noise, it is important to work with the write laser 
producing a Stokes photon with n < 1 (for example, n = 0.1 — 0.2). In Refs. [4]]], [3] an 
additional phase matching condition k„, + k s = k r + k as , where h w , k r are wave vectors of the 
write and read laser fields, and k s , k as are wave vectors of Stokes and anti-Stokes photons, 
respectively, is applied to make the axis of the write-read laser beams different from the axis 
of the scattered Stokes-anti-Stokes photons. Any Stokes photon that does not satisfy the 
condition of Raman resonance (or phase matching condition) and that, on the contrary, has 
the resonant frequency u es for the transition e — * s, will be spectrally broadened since it 
is produced due to spontaneous decay of the excited state |e) of a single particle, and not 
collectively. In our paper we show that such a photon or a photon produced by a single 
particle decaying in free space vacuum modes interacts differently with two- and three-level 
atoms. Therefore, the design of any quantum network should take into account the spectral 
properties of such a single photon and its interaction with atomic ensembles. 
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I K . 1: (a) Decay scheme of the excited particle emitting photons a{t) and b(t) in the cascade 
h — > e — ► g. (b) Detectors 1 and 2 detect photons a(t) and b(t), respectively. Detector 1 starts 
the clock when it detects a photon a(t) and detector 2 stops the clock when the photon b(t) is 
detected. The delay time between the two counts is stored and after many such events the detection 
probability of the photon b(t) versus time is reconstructed. 
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FIG. 2: Time dependence of the probability amplitude of a photon, which is formally presented 
as a sum of two parts. They are time domain counterparts of the symmetric and antisymmetric 
components of the photon spectrum (see the text). Plot (a) shows the time dependence of the 
symmetric component counterpart, and (b) - of the antisymmetric component. Dashed line: no 
absorber, but a photon is detected at distance I, so t is substituted by r = t — l/c. Solid line: with 
absorber of effective thickness T = 10. 
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FIG. 3: (a) Time dependence of the probability amplitude of the photon (dots) and its antisym- 
metric part (solid line) at the output of a thick absorber with optical thickness otol/T = 10. (b) 
Thickness dependence of the time integrated intensity of the symmetric (dash-dot line) and anti- 
symmetric (thick solid line) parts of the photon. Both are normalized to half of the time integrated 
intensity of the input photon C/q(0)/2. The thin solid line shows the limit of Beer's law exp(— 2Tb). 
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FIG. 4: Energy diagram of the absorber whose excited state e is coupled by to a metastable 
state m. Initially the particle is in the ground state g and a photon excites the transition g — > e. 
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FIG. 5: Time dependence of the function <fi+(l, r) for A p h/A e g equal to 0.1 (solid line) and (dots). 
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FIG. 6: Time evolution of the probability amplitude of the photon, b to t(l,T), transmitted through 
a sample of effective thickness = 30 (bold line). The thin solid line shows the same evolution 
with no absorber, 6o(0,r). The dashed line shows the probability amplitude, boi(l,T — td), which 
does not take into account the spectrum narrowing of the photon. The parameters of the absorber 
are V = 10j m and Q = 2T. The spectral halfwidth of the photon A p h is j m (a) and T (b). 
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FIG. 7: Time evolution of the probability amplitudes of the symmetric (a) and antisymmetric 
(b) parts of the photon without absorber, 6 S)(I (0,r) (dashed line), and with absorber of effective 
thickness T5 = 30, b s , a (l, r) (solid line). The parameters of the absorber and photon are the same 
as in Fig. 6(a). 
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